Gamow shell-model calculations of drip-line oxygen isotopes 
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We employ the Gamow shell model (GSM) to describe low-lying states of the oxygen isotopes ^'*0 
and ^''O. The many-body Schrodinger equation is solved starting from a two-body Hamiltonian de- 
fined by a renormalized low-momentum nucleon-nucleon (NN) interaction, and a spherical Berggren 
basis. The Berggren basis treats bound, resonant, and continuum states on an equal footing, and is 
therefore an appropriate representation of loosely bound and unbound nuclear states near threshold. 
We show that such a basis is necessary in order to obtain a detailed and correct description of the 
low-lying 1^ and 2^ excited states in '^^O. On the other hand, we find that a correct description 
of binding energy systematics of the ground states is driven by proper treatment and inclusion of 
many-body correlation effects. This is supported by the fact that we get ^^O unstable with respect 
to ^^O in both oscillator and Berggren representations starting from a ^^O core. Furthermore, we 
show that the structure of these loosely bound or unbound isotopes are strongly influenced by the 
^So component of the NN interaction. This has important consequences for our understanding of 
nuclear stability. 

PACS numbers: 21.10.-k, 21.30.-x, 21.60.-n, 24.30.Gd, 27.30.+t 
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Introduction. The study of nuclei far from stability 
is a leading direction in nuclear physics, experimentally 
and theoretically. It represents a considerable intellectual 
challenge to our understanding of the stability of matter 
itself, with potential implications for the synthesis of ele- 
ments. An important aspect of this research direction is 
to understand how magic numbers and shells appear and 
evolve with increasing numbers of neutrons or protons. 
Except for a qualitative understanding of the neutron 
and proton dependence of the magic numbers in terms 
of mean-field models, we lack a quantitative theoretical 
understanding in terms of the basic constituents of the 
underlying nuclear many-body Hamiltonian. 

For mesoscopic systems like nuclei, the interpretation 
of the standard magic numbers 2, 8, 20, 28, 50, 82, and 
126 is linked to the interplay between the filling up of 
single-particle (s.p.) orbitals (s.p. picture) and the un- 
derlying interactions. Magic numbers lead to so-called 
shell gaps in the s.p. spectra near the Fermi energy. For 
various magic numbers, in particular for stable doubly 
magic nuclei, one can have a large number of stable iso- 
topes (increasing number of neutrons TV but fixed number 
of protons Z) or isotones (increasing number of protons 
Z and fixed A^). Eventually, as one adds more neutrons 
or protons, viz., moving away from the valley of nuclear 
stability towards the drip lines, the outermost nucleons 
literally start to drip off the nuclei, thereby defining the 
very limits of stable matter. 

Nuclei near the drip line are open-quantum many- 
body systems for which the coupling with the scat- 
tering continuum plays an important role, and should 
therefore be explicitly taken into account. Configu- 
ration interaction (shell-model) methods, such as the 
Gamow shell model (GSM) (see Ref. [l| for a recent 



review of the GSM) or the continuum shell model 0], 
have been developed in order to properly include the 
coupling with the scattering continuum. The inclu- 
sion of continuum states complicates the solution of 
the many-body problem considerably, as the number of 
many-body basis states will explode. This has there- 
fore motivated different approaches to the description 
of loosely bound and unbound nuclear states, such as 
the Density-Matrix-Renormalization-Group (DMRG) Q 
and Coupled-Cluster approaches Furthermore, to 

complicate the solution strategies of the nuclear many- 
body problem is the need to include three-nucleon 
forces (3NFs). Recent results for light nuclei using ah- 
initio methods such as the Green's function Monte-Carlo 
(GFMC) and the no-core shell model (NCSM) [l], 
demonstrate that 3NFs are needed in order to give the 
correct binding energies and spectroscopy when compar- 
ing to experiment. It is still an open research problem 
what the role of 3NFs is in medium mass and neutron- 
rich nuclei close to the drip line. This applies also to our 
understanding of how different parts of the underlying 
nucleon-nucleon (NN) interactioii, such as the spin-orbit 
force 0, [1] and the tensor force 0, , affect the struc- 
ture of nuclei close to the drip line. 

In this work we study the ground- and low-lying states 
of ^^O and ^^O within the GSM framework. We derive 
for the first time a realistic effective shell-model interac- 
tion that includes the effect of the scattering continuum. 
The effective interaction is derived using many-body per- 
turbation theory (MBPT) starting with a Hamiltonian 
which reproduces NN scattering data. The choice of the 
oxygen isotopes is motivated by several reasons. First, 
the oxygen isotopes are the heaviest isotopes for which 
the drip line is well established. There are large exper- 
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imcntal campaigns worldwide (ill , [l^ . [l^ which aim at 
uncovering the properties of the oxygen isotopes, both at 
or close to the drip line and beyond. Two out of four 
stable even-even isotopes exhibit a doubly magic nature, 
namely ^^O {Z^8,N=U) [Hi and ^^Q {Z=8, N=16) 
[lE,[la|- The isotopes 25-28 q g^^,^ g^^l beheved to be unsta- 
ble towards neutron emission, even though ^^O is a dou- 
bly magic nucleus within the standard shcU-modcl pic- 
ture. Secondly, many-body descriptions starting with ef- 
fective two-body Hamiltonians defined for a model space 
comprising the lsi/2, 0^5/2, and s.p. orbitals, have 

always defied a p roper reproduction of the experimen- 
tal data [H, [3 HOl- To remedy this, the two-body 
shell-model interaction is usually fitted (fully or partially) 
0, [1^, [1^, im to reproduce experimental data for nuclei 
in the range 16 < A < 40. Without fitting the two-body 
Hamiltonians, which are typically derived using MBPT, 
shell-model calculations of the oxygen isotopes fare rea- 
sonably up to ^°0, but produce too compressed spectra 
in ^^O and lead to neutron separation energies that are 
always positive up to ^^O. All isotopes 1''-28q ^j^^^^ ^^g^j^ 
be reached via this model space result as well bound, in 
clear contradiction with experiment. These results per- 
tain to all available models of the NN interaction [3, H^l 
within the above model space. The ground states and 
excited states of the odd-even isotopes are also poorly 
reproduced (see for example Refs. [Tsl. ITqI). 

Recent ab-initio coupled-cluster calculations of the 
oxygen isotopes [l^l showed that chiral NN interactions 
can produce a rather flat binding energy curve of the oxy- 
gen isotopes ranging from to ^^O. It was illustrated 
that the existence of ^^O cannot be ruled out from ab- 
initio theory starting from modern chiral interactions. 
It was concluded that 3NFs will eventually decide the 
matter. Recent shell-model calculations of the oxygen 
isotopes within the s-d shell showed that inclusion of ef- 
fective 3NFs give added repulsion in the heavier oxygen 
isotopes and resulted in an unstable ^^0 [1^. 

Obviously, the understanding of the evolution of bind- 
ing energies and the location of the drip line in the oxy- 
gen isotopes are still unsolved theoretical problems. As 
shown in Refs. [2^. [23j. proper treatment of many-body 
effects are crucial in the oxygen isotopes. However, it 
is still an open issue what the role of the continuum is 
in the evolution of shell structure and in a correct de- 
scription of binding energy systematics in neutron-rich 
oxygen isotopes. As found in Ref. [11], a proper treat- 
ment of the continuum was necessary in order to explain 
the parity inversion of the ^^Be ground state. It is there- 
fore important to investigate what the role of continuum 
is on ground and excited states in neutron-rich oxygen 
isotopes. It is the aim of this article to throw light on, 
and disentangle the different effects coming from many- 
body correlations and the proximity of the continuum 
on the evolution of binding energies and shell structure 
in neutron-rich oxygen isotopes. Further, we will inves- 
tigate how different components of the NN interaction 
influence the low-lying resonance states in drip-line oxy- 



gen isotopes. 

Method, interaction, and model space. In this work we 
use the GSM to describe the low- lying states in ^''O and 
the ground state of ^^O. We construct the effective shell- 
model interaction starting from ^^O as a closed-shell core 
and use MBPT through second order The isotope 
^^O has been found to be a closed-shell nucleus with a 
considerable shell gap, making it suitable as a starting 
point for shell- model calculations. With ^^O as a closed- 
shell core we have at most three valence particles, making 
effects coming from a proper treatment of the continuum 
more transparent. So far, the effective interaction has 
typically been constructed using oscillator states as in- 
termediate particle states. In this work, we investigate 
for the first time what the role of the continuum is in con- 
structing a realistic effective shell-model interaction and 
in describing spectra of neutron-rich oxygen isotopes. 

The major advantage of the GSM is that it unifies 
structure and reaction properties of nuclei, that is, bound 
and scattering degrees of freedom are treated equally. 
Nuclei close to the scattering threshold have low-lying 
states that are in the continuum, and the GSM provides 
a simple framework for calculating the energy and life- 
times of these states. The starting point of the GSM is 
the Berggren completeness relation This s.p. basis 
treats bound, resonant, and scattering states on an equal 
footing, and it is the basic idea of the GSM to expand the 
many-body wave function in Slater determinants built 
from this basis. The representation of the Hamiltonian 
will no longer be Hermitian but rather complex symmet- 
ric, and a diagonalization will yield bound state spectra 
in addition to resonant state spectra. Hence, the energy 
and lifetime of the many-body resonant state comes out 
directly from the approach. 

In constructing the s.p. Berggren basis, we start from 
the one-body Hamiltonian for a spherical Woods-Saxon 
(WS) potential (see for example Ref. H^). The Berggren 
basis is obtained by diagonalizing the one-body Hamil- 
tonian in a plane wave basis defined on a deformed con- 
tour _L+ in the complex momentum plane (see Ref. pTj 
for details). This method is also known as the Contour- 
Deformation-Mcthod. There is freedom in choosing the 
contour L'^ provided the Schrodinger equation can be 
analytically continued in the complex momentum plane. 
Fully converged results should be independent of the 
choice of the contour, and in the following we used two 
different contours, and LJ, as given in Ref. [l^l, to 
check that our results are converged. The parameters of 
the WS potential are given by Vq = 55 MeV, Vso = 29 
McV, R = 3.1 fm, and d ~ 0.45 fm, where the shorthand 
Vso ~ 29 stands for the spin-orbit term. These parame- 
ters are chosen such that the energies of the s.p. orbitals 
vOdc,/2, vlsij2, and i^0c?3/2 are close to the experimen- 
tal values obtained by taking the binding energy differ- 
ences of ^^O-^^O and ^^O-^^O. We obtain, in units of 
MeV, e,od5/2=-6.53(-6.85), e,i,i/2 =-2.73(-2.74) and 
£iy(i3/2=l-43 — O.Oli (1.26 — 0.2j), where numbers in paren- 
thesis are experimental values. We refer hereafter to this 



basis as a WS Gamow basis. In our shell-model study 
of the oxygen isotopes ^"'O and ^^O we will use the WS 
Gamow basis for the neutron partial waves si/2, ^3/27 
and while for protons and all other partial waves, 

we use the oscillator representation. 

In constructing the effective shell-model interaction for 
neutron-rich oxygen isotopes, we start from the intrin- 
sic nuclear Hamiltonian H — t — ^coM + V , where t de- 
notes the kinetic energy, and tcoM is the kinetic energy 
of the center of mass. For the interaction V , we em- 
ploy the NN interaction by Entem and Machleidt [H, [s^l 
derived from chiral effective field theory (EFT) at next- 
to-next-to-next-to-leading order (N'^LO) and with a 500 
MeV cutoff. This interaction still has considerable high- 
momentum components, making MBPT non-convergent. 
In order to construct the effective shell-model interaction 
VcQ using MBPT, we use a low-momentum interaction 
l^=Mow-fc, obtained by integrating out those momentum 
components above a certain cutoff A [2^ . Thereafter we 
generate V in the WS Gamow basis using the procedure 
outlined in Ref. [23]. It is well known that this renor- 
malization induces forces of higher rank. In a completely 
renormalizcd theory, all obscrvablcs are independent on 
the cutoff A, and therefore the dependence of our results 
on the cutoff A tells us about the missing many-body 
physics. In this work we start from a ^^O core and limit 
ourselves to two-body effective shell-model interactions 
only. Since we have at most three neutrons outside the 
closed-shell core, it is reasonable to expect that the effect 
coming from effective three-body forces is small. 

In order to proceed with the construction of the effec- 
tive shell-model interaction, we need to define a model 
space for the valence neutrons. Since we will perform 
our GSM calculations starting with ^^O as a closed shell 
core, i.e., the 6^5/2 shell is filled, our many-body states 
are constructed by letting all valence neutrons act in the 
model space defined by the the s.p. orbitals Isi/2 a-nd 
ic?3/2, where i = 0, 1, . . . denote the i-th discrctizcd con- 
tinuum state for the partial wave. The effect of the 
continuum from the si/2 and ^5/2 partial waves are taken 
into account through the construction of the effective in- 
teraction. The effective shell-model interaction Vos is 
then constructed for this model space using MBPT [20| . 

Heff^Q-Q'j Q + Q' j qJq----. 

Here Q{uj) is given by diagrams which are valence-linked 
and irreducible. The argument lo is the starting energy, 
which is the sum of the unperturbed s.p. energies. The 
object Qiuj) is defined as 

Q{uj) = PHiP + PH,Q L^Qi/iP. 

In this work we take into account diagrams up through 
second order in perturbation theory in order to construct 
Q{uj). In Fig. [T]the diagrams that contribute up to sec- 
ond order in the effective interaction are shown. The 
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FIG. 1; The second-order contributions to Via. 



intermediate s.p. states in each diagram in can be either 
resonant, non-resonant continuum, or bound states. In 
our calculations we neglect the contributions from folded 
diagrams; see Ref. [1^ for further details. 

Results. We now turn to the GSM calculations of low- 
lying states in ^^O and ^^O. For the case of ^^O we have 
two neutrons outside the ^^O core, and we express the 
ground state of by filling the fsi/2 orbital. This 
orbital is well bound, and with a neutron separation en- 
ergy of S'„=2.7 MeV for [31|. We therefore expect 
mixing of higher-lying states to be small in the ground 
state of ^^O, and we take this effect into account through 
the construction of the effective interaction using many- 
body perturbation theory. For the excited states J = f 
and J = 2 of ^^O, we expand the wave function in a 
set of basis states |«J^) = |lsi/2 <8 id^/2]J^)- The or- 
bital denote the i-th discretized s.p. state of the 
WS Gamow basis, and z=f, • • -ngamow, where rigamow is 
typically between 20 and 35. In our calculations we used 
'T-gamow=30, and we checked that this was sufficient to 
reach convergence. 

In Fig. [2] the excitation energies of the low- lying states 
in ^''O and the ground state in ^^O with respect to the 
ground state of '^^O are shown. In these calculations we 
employed a cutoff A — 2.1fm^^ for the V|ow-fc interac- 
tion. We have varied the cutoff A around 2.1 fm~^ and 
found that our results are nearly independent of the cut- 
off A, meaning that the effect of neglected three- and 
many-body forces are small and will not change our con- 
clusions. Our GSM calculations give energies for the 
low-lying states that compare well with the experimental 
values. The splitting between the I"*" and the 2"*" states 
of ^^O compares well with experiment. The theoretical 
spacing obtained with our GSM approach is 540 keV, 
while the corresponding experimental value is 600 keV. 
Our calculated widths for the f + and 2+ resonant states 
are 0.44 and 0.26 MeV, respectively. In Ref. [H, Hoff- 
man reported that the widths are 0.03lg'j3 MeV for the 
1+ state and 0.05lo asMeV for the 2+ state. Our calcula- 
tions overestimate the width of the I"*" state as compared 
to experiment, on the other hand our calculated width 
of the 2"*" state is within the experimental uncertainties. 
As a comparison, we show results obtained with other 
effective interactions. The new universal sd Hamiltoni- 
ans (USDA and USDB), obtained by revising the USD 
interaction in order to be more suitable for exotic nu- 
clei jTot . results in a good agreement with experimental 
data. However, one should note that some effects coming 
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from the continuum and 3NFs are implicitly included in 
the new USD two-body matrix elements due to the fit 
made to experimental data. Although the USD interac- 
tion reproduces nicely two of the three experimental data 
points, the 2+ state in ^^O comes out as a bound excited 
state, a fact which disagrees with recent experiments. 
We also show results obtained with a pure harmonic os- 
cillator basis. In this case, the discretized continuum is 
not included in the model space, and we construct the 
effective interaction using perturbation theory through 
second order (see Fig. [T]), and a model space consist- 
ing of neutrons in the orbitals (lsi/2) 0^3/2)- Using a 
model space consisting of only bound oscillator orbitals 
puts the 1+ and 2+ states in ^^O too high in energy 
when comparing with experiment. With the inclusion of 
the continuum in the model space, the l"*" and 2+ states 
get lower in energy and also closer to experiment. This 
clearly shows an effect coming from a proper treatment 
of the continuum, and should therefore be included in a 
realistic description of low-lying states in ^^O. Further, 
the 2+ state in is a resonant state and cannot be 
properly described in an oscillator basis. 

Turning to the case of the ground state of ^^O, our 
GSM calculations correctly predict that it is particle un- 
stable, but puts it at a higher energy (£'=2.1, r=0.5 
MeV) as compared to recent experimental data (£=0.77, 
r=0.17 MeV). To get a better description of this state, 
configurations like i^(lsi/2)(«rf3/2)^ should be included 
since the GSM is currently solved in a filling configu- 
ration. Again we compare our GSM result for ^^O with 
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FIG. 2: The excitation energies of the J"" = 1+ and J"" = 2+ 
states in and the ground-state energy of ^^O measured 
from the ground states in ^'^O. 

a shell-model calculation where only the bound oscillator 
orbitals (lsi/2, 0^3/2) define the model space. An im- 
portant point to note is that ^''O comes out higher in 
energy as compared to ^^O, reflecting that it is unstable 
towards neutron emission. Even though the oscillator 



basis cannot provide us with a width or a lifetime for 
the ground state of ^^O, it predicts ^^O as particle un- 
stable and illustrates that the continuum coupling is not 
a dominant mechanism in explaining the position of the 
drip line of the oxygen isotopes. Thus, irrespective of 
whether we use a complex Gamow s.p. basis that defines 
a model space consisting of the (ls]^/2i ^^^3/2) states (with 
i ~ 0, . . . , rigamow), or a harmonic oscillator basis with the 
(lsi/2, 0^3/2) states only, ^^O is unstable against parti- 
cle emission as long as we use ^^O as a closed-shell core. 
This has important consequences for our theoretical un- 
derstanding of the stability of matter. As pointed out 
in a recent calculation by Hagen et al 23, many-body 
forces such as three-body interactions are crucial in ex- 
plaining the correct position of the drip lines. The fact 
that a calculation of ^^O with ^^O as a closed-shell core 
gives the correct qualitative description of the ground 
state of ^^O, hints at the fact that with three valence 
neutrons only, missing many-body physics does not play 
a central role. This should be contrasted to realistic cal- 
culations that utilize ^^O as a closed core. With nine 
valence neutrons, ^^O is strongly overbound using the 
standard model space comprising the oscillator orbitals 
(lsi/2, 0^5/2, 0^3/2) and a two-body effective interaction 
derived using MBPT. However, the location of the ex- 
cited states in ^^O and the width of the ground state of 
^^O depend crucially on the choice of basis, and we have 
shown that a GSM basis is needed in order to explain 
these states. 

We conclude this work with an investigation of which 
parts of the NN interaction are crucial for an under- 
standing of the above excited states in ■^"'O. To achieve 
this, we single out the ^Sq partial wave component of 
the NN interaction and vary its strength. This par- 
tial wave is particularly important to our understand- 
ing of pairing correlations in nuclei. In order to inves- 
tigate the importance of this partial wave on the low- 
lying resonance states in ^^O, we vary the strength as 
V{a)=a X (S'=0, £=0; J=0|y|5'=0, £=0; J=0), a e [0, 1] 
in a bare NN potential. Clearly, a=0 corresponds to 
an NN interaction with no contribution from this par- 
tial wave while a~l corresponds to the original interac- 
tion. The low-lying resonance states are calculated for 
each V{a), and the results are shown in Fig. [3l One 
can clearly see how the ^5*0 partial wave in the NN in- 
teraction influences the many-body resonances in ^"^O. 
The ground state gains additional binding when a in- 
creases since the ground state can be well described by 
the i^(lsi/2)^ configuration. At a=0 the position of l"*" 
and 2"*" are interchanged with respect to experiment, and 
they are degenerate if one analyzes the real part of the 
resonance energies. When a increases to the full interac- 
tion strength, the 1+-2+ splitting increases to 530 keV, a 
result close to experiment. The ^Sq partial wave, a cen- 
tral contribution to pairing in nuclei, is thus crucial in 
order to explain correctly the low-lying resonant states 
in 240. 

Conclusions. We have described the low-lying states in 
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FIG. 3; The energies of low-lying resonance states in ^''O, as a function of the strength a of the ^So partial wave component 
of the nucleon-nucleon interaction. See text for further details. 



the oxygen isotopes ^'^O and ^''O within the framework 
of the GSM and effective interactions derived through 
many-body perturbation theory. We have shown that 
there is a fine interplay between many-body effects and 
the proximity of the continuum. We demonstrated that 
the inclusion of a discretized continuum and resonant 
states are crucial in order to describe the excited states 
of ^^O and the ground state of ^'^O. A Gamow basis is 
needed in order to reproduce both the correct excitation 
energies and widths of these states. In order to predict 
the location of the drip line, we found that a proper inclu- 
sion of many-body forces and correlation effects are the 
deciding factors. Furthermore, the ^5*0 partial wave com- 
ponent of the NN interaction, an important contribution 



to strong pairing correlations in nuclei, plays an impor- 
tant role in explaining the correct value of these states. 
We conclude that the location of the drip line and excited 
states in the oxygen isotopes can only be determined by 
proper inclusion of both many-body correlations and the 
coupling with continuum degrees of freedom. 
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